Study of the ij{2S) decay to ppn^T: (tt^) 
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The branching ratios of ip{2S) —>■ pp-K^ir^ and ip{2S) — > ppuj, uj tt+tt^tt^ are measured 
and the tt'^'tt" invariant mass distribution in the first decay is discussed by analysing 14 miUion 
produced ip{2S) events collected by the BESII detector at the BEPC. 

PACS numbers: 13.20. Gv, 13.65.+i, 13.75.Lb, 14.40. Cs. 



^ ■ In the quarkonium model, the J/ip and ipi^S) are respectively the ground state'^1 and the first radial 
^ excitation of the cc bound statet^l As such, their decay is supposed to have a similar feature. In 
CN !both J/ip decay and ^ji^S) decay, the ppuj and ppp° are allowed three body states. For J/?/' decay, the 
Ch production of ppuj is observed at the fair fraction!^] of (1.30 ± 0.25) x 10~^ relative to other allowed 
1^ three body states, while only a upper limit is set to the branching ratio''^' of the production of ppp^. 

Therefore, to observe how the ppuj ppp^ are produced in ip{2S) decay is very interesting. 
^ . The CO decays to vr+vr" and 7r+7r~7r° at 1.7% and 89.7% branching ratios'^'^l The ppu is observed by 
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BESI^^l In the paper, we present the vr+vr mass spectrum in tp{2S) decay to ppn^n and the vr 
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mass spectrum in 'ip{2S) decay to ppvr+Tr 7r°. The analysis is based on 14 million produced ip{2S) events 
toUected by the BESII detector at the BEPC. 

This study shows an unknown bump of about 20 MeV width at 727 MeV in tt~^tt~' mass spectra. 
The BES apparatus is a magnetic spectrometer working at e~^e~ colliding mode, which has been fully 
described elsewhere ^^1. 

! The decay of ip{2S) to pp7i~^7T~TT^ is discussed fisrt beacuse of u decays to 7r"'"7r~7r° in large branching 
>< |ratio. 

Candidates for 'ip{2S) ppT^^-K^'K^ events are selected by requiring exactly four reconstructed 
charged tracks in the drift chamber with zero net charge. Tracks with transverse momentum p^y > 0.08 
,GeV and | cosOch \< 0.85 are accepted, where 6ch is the polar angle with respect to the beam direction. 
The particles are identified by requiring that their combination weights of time-of-fiight (TOP) and 
the ionization energy loss (dE/dx) in the drift chamber be consistent with the corresponding particle 
hypothesis. The events with at least two charged particle satisfying proton and anti-proton hypotheses 
■and one satisfying pion hypothesis are selected. 

To remove the contamination from ip{2S) —* AAtt^, Mp^- > 1.15 GeV and Mp^+ > 1.15 GeV are 
required. 

An isolated photon is defined as a cluster in the barrel shower counter with at least two readout 
layers hit, energy larger than 30 MeV, outside a 25° cone around the p and outside a 12° cone around 
other charged particles to reject possible fake photons produced by p annihilation and/or radiated 
by other charged particles inside the shower counter. The cluster is also required to has its incident 
direction (from the interaction point to the first hit point in BSC) inside a 20° cone around the cluster 
developing direction. 

The events are kinematically fitted for the ip{2S) ppTi^-K^'K^ topology by imposing energy and 
momentum constraints (4C). Xpp^+^-^T, < 20 and < to remove ^(25") ppn+n-'y. 

The combination with the smallest the 4C fit is chosen to identify the radiative photon if there 
are more than two photon candidates in an event. Fig. 1 showes the 77 invariant mass distribution. 
The 7r° signal is seen at 0.135 GeV. To select 7?°, it is required | M^^ - 0.135 |< 0.040GeV (3a). 

A contamination is 4'i2S) — > 7i~^7i~J/ip, J/ip ^ ppn^. To remove it it is required | M^p^o — 3.097 |> 
0.065 GeV (3cr). (why do not use Mpp^o < 3.0 GeV, because Mpp^o distribution is around 3.0 GeV for 
ip{2S) ppuj, UJ TT^Tr^Tr".) 
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Fig. 1. 77 invariant mass for ip{2S) — > ppTT+vr 77. 



550 events are obtained for ip{2S) — > ppTT~^n 7r°. tt+tt tt^ invariant mass is shown in Fig. 2. At 
0.783 GeV is u signal. The mass resolution and efficiency though Monte Carlo simulation are obtained 
to be 14.6 MeV and 3.6%. 
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Fig. 2 7r"'"7r 7r° invariant mass for ip{2S) — >• pp tt+tt 7r° 



The a; mass, width, event number can be obtained using unbinned Breit-wigner fit. The background 
is fitted with 4th order nomial. = 783.4 ± 4.5 ± 0.2 MeV. = 12.7 ± 7.1 ± 3.0 MeV. Nevents = 
35.0 ±6.0 ±6.7. 

There are not contriburion to u signal from continuum data, which id derived from 6.42(1 ±4%)p6~^ 
data at 3.65 GeV. 
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The branching ratio is Br{ip{2S) — > ppu,uj tt+tt 7r°) 
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= (6.9 ± 1.2 ± 1.3) X 10- 



Br{ip{2S) ppu) = (7.7 ± 1.3 ± 1.4) x 10"^ 

The BES I data gives'^l Br{ij{2S) ppu) = (8.0 ± 3.0 ± 1.0) x 10"^ 

Now the decay of '4>{2S) to ppn'^n'' is discussed. The charged particle are selected as above. To 
remove the one or multiple photon events, the missing momentum is required to be Pmissing < 0.1 GeV. 
■0(25') A A is a contamination source. It is removed by Mp^4 



> 1.15 GeV and Mp^- > 1.15 GeV. 



Another contamination source is 'ijj{2S) tt+tt J/i/j, J/ip pp. It is removed by requiring Mpp < 
The ppn^n^ mass is shown in Fig. 3. It is '4'{2S) signal. Its resolution is 36 MeV. 
- 3.686 |< 0.1 GeV is required. 



3.0 GeV. 
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To remove ip{2S) ppK+K 4c fit is performed. xl{2S)^pp7r+n- < 30 and x|(2S)-^-~+~- < 
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X^(2S)^ppK+K- is required. 
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Fig. 3. ppTr'^Tr~ invariant mass distribution for ip{2S) ppK'^n' candidates. 

The TT+TT^ invariant mass distribution is shown in Fig. 4. 1691 events are there. Two narrow 
bumps can be seen at 727.0 and 783.0 MeV and a httle bump at 977.0 GeV, i.e. the position of /o(980). 
From Monte Carlo simulation, the BES-11 detector acceptance efficiencies and the mass resolutions at 
three points are obtained. The mass resolutions are 6.0 MeV at the first two points and 6.8 MeV at 
977.0 GeV. The efficiencies are 21.4%, 21.8% and 22.4%. 

Their masses, widths, event numbers and statistic significance (S.S.) can be obtained using unbinned 
Breit-wigner fit. The background is fitted with 3th order nomial. Each of the first two bumps is described 
using free parameters: area, width, mass and one fixed parameter: resolution. The mass and width of 
/o(980) are fixed at the values of PDG2004 977.0 MeV and 44.0 MeV. The mass resolution of /o(980) 
is fixed. 

The table 1 shows these values 




Table 1. The parameters of bumps by fitting the first two as unkown ones. 





mass (MeV) 


width (MeV) 


Event 


S.S. (a) 


bumpl 


726.9 ±4.1 ±2.1 


20.7 ± 12.5 ±8.1 


69.0 ± 15.8 ± 12.8 


4.9 


bLinip2 


783.0 ±2.2 ±0.1 


3.2 ± 5.6 ± 1.2 


12.5 ± 10.6 ± 10.7 


5.0 


/o(980) 


977.0 


44.0 


8.8 ± 16.41^;^ 


0.5 



statistic significance is defined as significance = y2 x {InLmax — InL^) Where L^ax is hkehhood 
function value for 1. the first two signals yield with resolutions fixed and areas, masses and widths 
to float, 2. the /o(980) signal yield with resolution, mass, width flxed and area to float and 3. the 
background shape parameters to float. LO is the likehhood function value for the signal hypothesis 
under consideration corresponding to a zero yield and the other two signals yield as in Lmax- 

The mass and width of bump2 are consistent with those of u. Bumpl is unknown. The known 
particles listed in PDG2004W have no consistent masses and widths with those of bumpl. 

The branching ratios are 

Br{ip{2S) ppbumpl, bumpl n'^ir') — 
^'"""'"•"^ = (2.3 ± 0.5 ± 0.4) X 10-^ 

Br{ip{2S) ppbump2, bump2 tt'^tt") 
= (1.4 ±0.4 ±0.4) X 10"^ 

5r(^(25) ^pp/o(980),/o(980) ^ tt+tt") 
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= (2.8 ±5.2+^1) X 10"^ 

From Br{ip{2S) ppuj,uj 7i^7T~tt^), it is derived Br{ilj{2S) ppu!,uj tt+tt^) = (1.3 ± 0.2 ± 
0.2) X 10^^, which is not consistent with branching ratio of bump2 production. 



Br{i){2S) PP71+71-) = ^ r.^n., 6 = (^-O ± 0.1 ± 0.6) X 10" 



(1691 - 123 - 27) ± 40.0 ± 175.6 
0.22 X 14 X 10^ 

where AA is excluded, 123 is contriburion of continuum data, which id derived from 6.42(1 ± 4%)p6~^ 
data at 3.65 GeV. 27 is from sidebands. 40.0 is statistical error. 175.6 is systematic error. The efficiency 
is 22.0%. The sideband events are obtained using the cut 0.1 <| MppT^+j,- — 3.686 |< 0.2 GeV and the 
same other cuts for Fig. 4. 

In addition to AA, two decay processes of ip{2S) have ppn^n^ final state. One is ipi^S) A++A 
with the branching ratio of (1.28 ± 0.35) x 10~^ given by Another is direct decay ip{2S) 

ppn+n-. The branching ratio given by PDG2004[4'9l is Br{^lj{2S) ppn+n-) = (8.0 ± 2.0) x 10-^ 
which includes the contribution from the AA, and direct decay. 

Then what is branching ratio of direct decay ip{2S) — > ppTT^Ti~ from this measurement? Using 
Br{%l){2S) —>■ A^+A ), Monte Carlo simulation shows that 368 A^+A events contribute Fig. 4, i.e. 
number of direct decay ppir'^Tc' is 1691-123-27-368-65.9-42.5=1064.6, two bumps are excluded (because 
number of /o(980) is small and has large statistics, it is ignored). The corresponding branching ratio is 
3.5 X 10~^. In Fig. 4 the solid line is from Monte Carlo evevts of ip{2S) A^^A with branching 
ratio of 1.28 x 10"'^ plus direct decay ip{2S) — > ppTT^Ti~ with branching ratio of 3.5 x 10~^. 
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Fig. 4. vr+vr" invariant mass distribution for ip{2S) — > ppvr+vr". The shaded area is data. The 
smooth line is from BW fit. The solid line is from Monte Carlo simulation ip{2S) — > ppTi^Ti~ and 
ip{2S) A"'""''A . The small peak at 0.782 GeV near bottom line is from Monte Carlo simulation 
with Br{ilj{2S) ppu,uj ^ tt+tt~) = 1.3 x 10"^ 

The sample of 14 x 10^ Monte Carlo inclusive il'{2S) decay from lund-charm generatort^^l is used 
to estimate the background. The sample is throught the same cuts. 4686 events are selected. Among 
the events, the contamination is shown in table 2. Monte Carlo simultion for each channel shows these 
channels maybe have no contributions to ppvr+vr" invaraint mass distribution around 0.727 and 0.783 
GeV. 



Table 2. Contamination from 'ip{2S) decay to anything 
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source 


number 


bran, ratio 


ti'tW, n' — > ^0 

' 1 J/J/1 ' 1 1 r 


22 




TIVV, T) — > 0'7r~'~7r~ 


13 




/vcl 1 7 A,cl Irirr 


1 


< 4.5 X 10"^* 




5 


< 9.5 X 10~^* 


yr27, yr2 ^ al(1260)°7r° 


1 






1 


9.5 X 10"^ 


7r°7r+pA — 


1 






1 


1.4 X 10""^ 




2 






1 


< 1.8 X 10^2 


sum 


48 





here, branching ratios come from PDG2004. *Br{xci1, Xci PPP) < Br{xci1,Xci PW^t^ ) — 
4.5 X 10-^ Br{xc2l,Xc2^PPp) < 5r(xc27, Xc2 ppvT+Tr") = 9.5 x 10"^ 

Systematic errors for eficiency are caused by difference between data and MC simulation. Studies 
have determined these errors to be 8% for the tracking cfficcincy, 2% for photon identification, t^*'', 5.6% 
for PID, 1.9%, 11.1%, 16.3%, —78.6%, 8.4% for kinematic fit oi ppn^Ti^ , ppbumpl, ppbump2, pp/o(980) 
andppa;, 0%, 0.01%, 4.1%, 62.5% and 1.2% for BW fit of ppn~^TT~ , ppbumpl, ppbump2, pp/o(980) using 
4th order nomial to describe the background and ppcj using 5th order nomial to decribe the background, 
1.1%, 1.7%, 8.5%, -12.9% and 5.9% for | cos9 \< 0.8 for ppn+n~, ppbumpl, ppbump2, pp/o(980) and 
ppu, 0.5%, 0.5% and 3.1% for cffcicncics of bumpl, bump2 and /o, 0.6%, 9.7%, 10.5% and -100.0% 
for Pmissing foT ppTT+TT", ppbumpl, ppbump2 and pp/o(980), 5.0% for number of ip{2S) events. Total 
systematic errors are 11.2%, 18.5%, 25.2%, IiHq and 15.4% for ppTr+Tr", bumpl, bump2, /o(980) and lu 
respectively. 

The parameters of the two bumps in Fig. 4 are obtained by fitting them as completely unkown 
partilces. But theoretically ^{2S) can decay to {J^'" = 1 ) and ppu (1 ). This gives rise to 
the questions: l.the and u plus their interference can fit? 2. the u plus one unknown bumpl can 
fit? 3. the p° and lv and their interference plus one unknown bumpl can fit the two bumps? 

Fig. 5 showes Breit- Winger fits for these three combinations, in which the the masses and widthes 
of the p^, OJ and /o(980) are fixed at PDG2004. The fitting parameters are given in table 3. 



Table 3. Breit- Winger fitting parameters 



combination 


bump 


mass (MeV) 


width (MeV) 


event 


S.S. 


Br.(10-^) 


1 




776.0 


149.0 


156.6 ±31.9 ±36.8 


5.1 


5.4 ± 1.1 ± 1.3 




782.5 


8.5 


13.0 ±7.7 ±5.1 


1.8 


0.43 ±0.25 ±0.17 


/o(980) 


977.0 


44.0 


20.8 ± 16.0 ±9.4 


1.3 


0.66 ±0.52 ±0.30 


2 








49.9 ±12.2 ±10.1 


5.0 


1.6 ±0.4 ±0.3 


/o(980) 






9.5 ± 16.4+^:^ 


0.6 


0.30 ± 0.52+^;^^ 


bumpl 


726.7 ±4.0 ±2.5 


19.4 ±10.8 ±7.0 


66.5 ± 15.6 ± 14.9 


4.9 


2.2 ±0.5 ±0.5 


3 








18.5 ±13.2 ±11.2 


4.5 


6.4 ±4.6 ±3.9 








35.8 ± 15.8 ±6.2 


3.8 


1.2 ±0.5 ±0.2 


/o(980) 






20.1 ± 15.6 ±9.7 


1.0 


0.64 ±0.50 ±0.31 


bumpl 


726.8 ± 4.2 ± 2.3 


18.4 ±12.5 ±7.2 


63.7 ± 14.4 ± 11.4 


4.0 


2.1 ±0.5 ±0.4 



5 






(1).M^+^- {GeV) (2).M^+^- {GeV) (3).M^+^- {GeV) 

Fig. 5 Three kind of combinations of p° uj /o(980) and bumpl to fit the tt+tt" invariant mass distribu- 
tion from ^'(25') — » ppvr+vr". (l)p° + uJ +p^io interference + /o(980). {2)unknown bumpl + u + /o(980). 
{3)unknown bumpl + p° + a; + pa; interference + /o(980). 

In this work, br{ip{2S) ^ ppn+n-) = (5.0±0.1 ±0.6) x 10"^ and 6r(^(25) ^ ppu) = (7.7± 1.3 ± 1.4) x 
10~^ are measured by analysis of BES-II 14 milhon ip{2S) events sample. This leads to the observation 
of two bumps in vr+vr" mass spectrum in the decay process '?/'(2S') — > ppvT+vr". The bumps are not 
observed in J/ip pp7r'^7r~W. The mass and width of the first bump can not match those of any 
particles in PDG2004. Another bump has a consistent mass and width with those of uj. But if bump2 
is real u, there is no consistency between the branching ratio for ip{2S) — > ppbump2, bump2 tt+tt^ 
and the branching ratio for ip{2S) — > ppu, cu Tr+Tr'vr". 

Duo to small statistics, the quantum number J^'^ is not set for bumpl and bump2. The bumpl has 
a statistical significance 4.0 < S.S. < 4.9 a. Whether the two bumps come from statistical fluctuation 
should be confirmed by larger statistics of 4'{2S) events. 
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